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ABSTRACT 

We present an observational reconstruction of the radial water vapor content near the surface of 
the TW Hya transitional protoplanetary disk, and report the first localization of the snow line during 
this phase of disk evolution. The observations are comprised of Spitzer-IRS, Herschel-PACS, and 
Herschel-HIFI archival spectra. The abundance structure is retrieved by fitting a two-dimensional 
disk model to the available star+disk photometry and all observed H 2 lines, using a simple step- 
function parameterization of the water vapor content near the disk surface. We find that water vapor 
is abundant (~ 10 -4 per H2) in a narrow ring, located at the disk transition radius some 4 AU from 
the central star, but drops rapidly by several orders of magnitude beyond 4.2 AU over a scale length 
of no more than 0.5 AU. The inner disk (0.5-4 AU) is also dry, with an upper limit on the vertically 
averaged water abundance of 10~ 6 per H2. The water vapor peak occurs at a radius significantly more 
distant than that expected for a passive continuous disk around a 0.6 M Q star, representing a volatile 
distribution in the TW Hya disk that bears strong similarities to that of the solar system. This is 
observational evidence for a snow line that moves outward with time in passive disks, with a dry inner 
disk that results either from gas giant formation or gas dissipation and a significant ice reservoir at 
large radii. The amount of water present near the snow line is sufficient to potentially catalyze the 
(further) formation of planetesimals and planets at distances beyond a few AU. 

Subject headings: planetary systems: protoplanetary disks — astrochemistry — stars: individual (TW 
Hya) 



1. INTRODUCTION 

Giant planets and planetesimals form in the chemi- 
cally and dynamically active enviro nments in the inner 
zone of gas-rich protoplanetary disks (jPollack et al.lll996t 
IArmitagal2011l ). Processes taking place during this crit- 
ical stage in planetary evolution determine many of the 
parameters of the final planetary system, including the 
mass distribution of giant planets, as well as the chemical 
composition of terrestrial planets and moons. 

Indeed, some of the most important features of mature 
planetary systems may be dictated by the behavior of the 
water content of protoplanetary disks. Such processes in- 
clude the formation of a strong radial dependen ce of the 
ice ab undance in planetesimals - the snow line (|Hayashil 
1 198 If ) - The deficit of condcnsible water inside the snow 
line leads to a need for radial mixing of icy bodies during 
the later gas-poor stages in the evolution of the plan- 
etary system, in pa rticular to explain th e existence of 
water on the Earth ([Raymond et alJl2004[ ). Ice mantles 
on dust beyond the snow line allow grains to stick at 
higher collisional velocities allowing for efficient coagu- 
lation well beyond the limit s of refractory grain growth 
alone ([Blum fe Wur m 2008). The combined mass of the 
solid component in protoplanetary disks is likely dom- 
inated by volatile, but condensible, species, the most 
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abundant of which is water ([Lecar et al.ll2006l ). 

Thus, core accretion models that rely on the avail- 
ability of large concentrations of solid material sug- 
gest that giant planets will initially be found beyond 
the snow line, as their forma tion requires the existence 
of a massive ice reservoir (Kennedy fc Kenvonl 120081 : 
iDodson-Robinson et al.l l2009f) . The added solid den- 
sity offered by condensible water may also offer a way 
out of the "meter-bar rier" for planetesimal formation 
([Weidenschillingl Il997| ). in which centimeter to meter 
size icy bodies migrate inwards and are lost to the cen- 
tral star on time scales mu ch shorter than tho se of 
steady state grain aggregation dBirnstiel et al.l2009l ): Be- 
yond a certain mass density of solids, local concentra- 
tions of boulder-sized particles may overcome the effects 
of turbulence acting to disperse them and will grav- 
itationally contract to form Ceres-mass planetesimals 
([Johansen et al.|[2007l [2009]) . 

Given the central role that water plays in the formation 
of planetary systems, its distribution and abundance in 
protoplanetary disk s has been the subject of intense the- 
oretic al study (e.g., iCiesla fc Cu zzi 2006; Ga raud fc Linl 
[2007tlKretke fc Linil2007D . Yet, the only source of obser- 
vational constraints on the initial radial distribution of 
water in planet-forming regions ha s been the dist ribution 
of ice in the current Solar System (|Havashilll981l) , the re- 
maining bodies of which represent only a small part of 
the story. Apart from the fact that the solar system wa- 
ter distribution may be unique, theoretical study shows 
that the initial distribution of ices in protoplanetary disks 
evolves as does the disk temperature distribution (e.g. 
iGaraud fc Lin|[2007t ). and that in later stages it is con- 
founded by the dynamical interactions between planets 
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and the planetesimal swarm (J Gomes "eFaIll2005h . The fi- 
nal distribution of water in a planetary system is thus an 
obscured tracer of initial conditions; and to understand 
the role of water and other volatiles in the formation 
of planets, we should measure their distribution before, 
during and after planet (csimal) formation. 

The most direct way to trace the water abundance dis- 
tribution in protoplanetary disks is to image line emis- 
sion from water vapor. However, due to the strong radial 
gradient of gas temperature in disks, each line of a given 
excitation energy will only trace a limited areal extent of 
the disk. For instance, high- J rovibrational transitions in 
the mid-infrared waveband are sensitive to the warmer 
(> 300 K) and inner parts of the disk surface (< 10 AU). 
Far-infrared transitions arise from cooler gas (< 300 K) 
in the outer (> 10 AU) and deeperregions of the disk. In 
order to measure the total water vapor content in the 
disk surface, it is necessary to observe lines across states 
of widely varying excitation energies. 

That line emission from water vapor is common in 
disks is now well established thanks to observations with 
the Spitzer IRS instrument, which has d etected a forest 
of em ission lines in the 10-35 /jm range. iCarr fc Na]ital 
(|2008D and lSalvk etaLl (|2008f ) first reported a large num- 
ber of rotational lines due to warm water in AA Tau, AS 
205N, and DR Tau, while iPontoppidan et all (|2010b( l re- 
port a water emission detection rate of ~50% in their 
Spitzer survey of 46 disks surrounding late-type (M-G) 
stars. With upper energy levels of E up ~ 1000 — 3000 K, 
these water lines arise from the inner few AU region of the 
disks. Detailed radiative transfer models show that the 
surface water vapor abundance may need to be truncated 
beyond ~1 AU for a typical classical T Tauri star (cTTs) 
disk in order to match the g eneral line ratios over 10- 
35 /urn (jMeijerink et al.l l2009t ). Follow -up ground-based 
spectra of the brightest targets have constrained the gas 
kinematics, and con firmed that the water vap or resides 
inside the snow- line ([Pontoppidan et aT1l2010al) . 

Recently, complementary spectral tracers of water va- 
por in the outer disk have become accessible via Her- 
schel Space Observatory PACS and HIFI observations. 
iRiviere-Marichalar et al.1 (|2012f ) describe the discovery of 
warm water emission at 63.3 /tm in 8 out of their 68 
T Tauri disk sources, and tentative detections of wa- 
ter transitions with similar excitation energies have also 
been reported for the Herbig Ae/Be star HD 163296 
(|Meeus et all 12011 iFedele et al J 120121 ). The first sen- 
sitive search for ground-state emission lines of cold wa- 
ter vapor in DM Tau indicated that the vertically av- 
eraged water vapor abundance in the outer disk of this 
source is extremely low, < 10~ 10 per hydrogen molecule 
(jBergin et al.l 120101) . Finally, two ground state water 
emission lines in TW Hya were detected with HIFI 
(|Hogerheijde et ahl 120111 ). Again, the vertically averaged 
water vapor abundance is < 10 -10 per hydrogen, while 
the estimated peak water vapor abundance is closer to 
^10~ 8 — 10~ 7 in the near surface layers of the disk. These 
are values that can be produced by UV photodesorption 
from icy dust grains. In this interpretation, the low ex- 
citation water lines trace a large, but otherwise unseen, 
reservoir of ice in the outer disk. 

Here we present an observationally-based method for 
reconstructing the water vapor column density profile 



in the surface of a protoplanetary disk, based on multi- 
wavelength, multi-instrument IR/submillimcter spectra, 
and apply the method to the TW Hya transitional disk. 
We chose TW Hya because of the availability of deep 
archival Spitzer and Herschel water spectroscopy and 
the existence of detailed structural models. Further, we 
report on the spectroscopic identification of warm water 
vapor emission at 20-35 /mi, a detection that permits, in 
combination with the Herschel data, the first localization 
of the snow line in a transitional disk. 

The outline of the paper is as follows: In Section [21 
we describe the Spitzer and Herschel data. In |j3l we 
describe the structural gas/dust model for the TW Hya 
disk. Section [4] applies a two-dimensional line radiative 
transfer model to the model disk structure to retrieve a 
radial abundance water profile based on the full spectro- 
scopic dataset. The results are discussed in Sj5] 

2. DATA REDUCTION 

The Spitzer IRS 10—35 /tin high resolution mode 
(R~600) spectra of TW Hya were acquired as part of 
a survey program of transitional disks (PID 30300), with 
J. Najita as PI. The observations were done in two 
epochs using different backgro und observat i on st rate- 
gies, for an observation log see iNaiita et al.l (f2010h . In 
the first epoch, TW Hya was observed on source (AOR 
18017792), followed by north and south off-set observa- 
tions (AORs 18018048, 18018304). In the second epoch 
(AOR 24402944), TW Hya was observed using a fixed 
cluster-offsets mode, such that the background scans 
were observed in the same sequence. All of the datasets 
were extracted from the Spitzer archive and reduced us- 
ing the Ca ltech High-resolutio n IRS P ipeline (CHIP) de- 
scribed in IPontoppidan et al.1 (j2010bl ). which takes full 
advantage of the existence of redundant background ob- 
servations. CHIP implemen ts a data reduction s cheme 
similar to that developed bv lCarr fc Naiital (I2008D . 

In an ex tensive analysis of t he 10-20 /mi Short-Hi (SH) 
spectrum. INaiita et al.l (|2010l ) report the detection of ex- 
cited OH emission lines and a host of other molecular 
features from TW Hya, but no water emis sion above a 
lower flux limit of ~10mJy. As stressed bv lNajita et al.l 
(|2010l ). the lack of detectable water emission from the 
highly excited water lines below 20 /tm does not preclude 
the presence of cooler water vapor, and in a follow-up pa- 
per (in prep), these authors interpret the SH and Long-Hi 
(LH) data together. Our independent reduction of the 
deep Spitzer IRS observations at longer wavelengths dis- 
plays numerous water emission lines in the 20—35 /tm 
LH module. The strongest transition at 33 /tm can be 
clearly seen in all LH epochs. The spectrum from the 
fixed cluster-offsets AOR has the highest signal-to-noise 
ratio (SNR), and our analysis is therefore based on this 
spectrum (presented in Figure [TJ). The SH spectrum is 
scaled by a factor of 1.25 to match the flux of LH spec- 
trum and IRAS 12 /tm photometry data. 

The Herschel HIFI line fluxes for the ground state 
ortho and para water transitions are taken from 
iHogerheijde et aLl (120111 ). Somewhat higher excitation 
lines are probed by the PACS instrument, and narrow- 
range line spectra of TW Hya from 63—180 /tm were ac- 
quired as part of the public Herschel Science Demonstra- 
tion Phase program (observation ID 1342187238). The 
PACS data were reduced using the Herschel interactive 
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processing environment (HIPE v. 6.0), up to level 2. After 
level 1 processing, spectra for the two nod positions were 
extracted separately from the central spaxel. The spectra 
of the two nod positions were uniformly rebinned, using 
an over-sampling factor of 2 and an up-sampling factor 
of 2, before co-adding to produce the final result. Linear 
baselines were fitted to the local continua and integrated 
line fluxes (see Table [2]) were obtained using Gaussian 
fits. While no water emission is formally detected with 
PACS at 3cr, there are indications of water lines consis- 
tent with our water abundance model predictions; see 
Section |4] for further discussion. 

3. THE TW HYA DISK STRUCTURE 

The accuracy of retrieved chemical abundances in pro- 
toplanetary disks is coupled to how well the two- or 
three-dimensional gas kinetic temperature and volume 
density structures are known. Our model development 
consists of three stages: First, we construct an axisym- 
metric model of the dust component of the disk, based on 
fits to the continuum spec tral energy distribution (SEP) , 
using the RADMC code (|Dullemond fc Pominikl 120041 ). 
We then assume that the gas density distribution fol- 
lows that of the small dust grain population — a rea- 
sonable assumption since the small grains carry most of 
the disk opacity and are dynamically well-coupled to the 
disk gas. However, we do allow thermal decoupling of the 
gas from the dust near the disk surface, in the form of 
a constant temperature offset constrained using the CO 
rovibrational emission spectrum near 4.7 /mi. Finally, we 
constrain the water vapor radial column density distribu- 
tion (and thus abundance) at all disk radii via fits to the 
multi-wavelength water spectra acquired for TW Hya. 

3.1. The dust temperature and density structure 

Dust, as the dominant source of disk continuum opac- 
ity, determines th e disk SEP shap e. At a distance of 
only - 51 ± 4pc (|Mamaiekl [2005f) . TW Hya is one of 
the best studied protoplanctary disks, with extensive 
photometric data a vailable from UV to cm wavelengths. 
ICalvet et al.1 ()2002h developed a physically self-consistent 
dust structure model for TW Hya. They found that 
the disk is vertically optically thin within 4AU, corre- 
sponding to a depletion of small dust grains by orders 
of magnitude at these radii. This inner zone of dust de- 
pletion has been confirmed by subseq uent observations 
(|Eisner et afll2006t iHughes et al.ll2007t ). A small amount 
of dust in the inner disk is needed to explain the 10 - 
25/Ltm amorphous silicate features (jCalvet et al.l 120021 ). 
and accretion onto the star continues. The existence of 
hot gas in the innermost disk has been further confirmed 
by the dete ction of CO Ay=l rovibrational emission 
near 4.7 /im (jRettig et alJl200l [Sarvk et al.ll2007l 120091) 
and via Sp ectro-Astrometric (SA) o bservations of these 
same lines ([Pontoppidan et all 120 08). It is not yet clear 
whether the gas content has been depleted to the same 
degree as the dust in the inner disk, although lGorti et al.l 
(|2011| ) estimated enhancements of 5-50 in gas/dust in in 
the inner disk. 

We adopt M* = 0.6 M Q , R = 1 M Q and T cff = 4000 K 
for the stellar mass, stel lar radius and effective temper- 
ature (fWebb et al.l ll999). The inclination ang le of the 
disk axis to the line-of-sight is fixed to i = 7° (|Qi et al.l 



Table 1 

Disk parameters for the dust model 



Parameter 


Inner disk 


Cavity wall 


Outer disk 


r (AU) 


0.055 - 4 


4 -5 


5-196 


P 


0.8 


0.0 


1.0 


7 


0.3 


-0.25 


0.2 



Note. — Here 7 is the index that describes the 
vertical scale height Zd- Zd = zo(r /ro) 1 ^ 7 where 20 = 
0.7 at 7*0 = 10 AU; p is the index of the surface density 
distribution: = Eo(r/ro)~ p . 



120041 ). The standard gas-to-dust ratio of 100 is used to 
estimate the total disk mass. 

The dust s t ructur e model adopted is similar to 
lAndrews et al.l (|2012t) : an optically thin inner disk be- 
tween the sublimation radius (r su b) and the cavity edge 
(r cav ), a cavity "wall" and an optically thick disk. The 
cavity "wall" component is needed to repro duce the spec- 
tral shape of TW Hya from 10—30 /im (jUchida et al.1 
120041 ). The sublimation radius was set to r su b = 
0.055 AU, the location where dust temperature reaches 
1400 K. The cavity radius is fixed as r cav = 4 AU, in ac- 
cord with previous near-IR interferometrv ([Eisner et al.l 
l2006HHughes et al.ll2007l : lAkeson et al.ll201lD. We ad opt 
the power-law density profile of lAndrews et al.l (|2012l ) for 
the outer disk, which is based on 870 /im continuum in- 
terfcrometric imaging. A summary of the dust structure 
model para meters can be s e en in Table L 

Foll owing iPollack et all (|1994l ) and iP'Alessio et~aT1 
(|2001h . we use a dust mixture of astronomical silicates, 
organics and water ice. For the optically thin inner disk, 
pure glassy silicate is used to match the strong silicate 
bands at 10 and 18 /im. The dust size distribution is 
taken to be n{a) oc a ~ s with s = 3 . 5 bet ween a m - ln = 0.9 
and a max = 2.0 /im (jCalvet et al.H2002T ). In the cavity 
wall, we use a dust size range between 0.005 and 1 /im; 
outside the cavity wall, 95% (by mass) of the dust has a 
size distribution extending from 0.005 /im to 1 mm, with 
the remaining 5% contained in grains from 0.005 - 1 /im. 
Total dust masses in the optically thin and thick parts 
of the disk are 2.77 x 10~ 9 and 4.4 x 10" 4 M . The SEP 
fit and our fiducial gas temperature and density profiles, 
discussed further below, are presented in Figure [2l 

3.2. The gas density and temperature 

The gas/dust ratio in disks evolves due to g rain growth 
and transport processes (|Birnstiel et al.ll2009l ). and likely 
has significant radial and vertical structure. For simplic- 
ity, the gas density is here assumed to follow that of the 
dust, with a constant gas/dust mass ratio throughout 
the disk. For TW Hya, th e estimated global gas /dust 
ratio varies from 2.6 to 100 (|Thi et al.ll2"010l : iGorti et al.1 
I2011h . The recent detection of the HP J = 1-0 line from 
TW Hya has offered an independent an d robust estima- 
tion of the gas mass (|Bergin et al.ll2013l ) , which indicates 
a globally averaged gas/dust ratio close to 100, a value 
we adopt here. Significant dust vertical settling will pro- 
duce an enhanced gas-to-dust ratio in the upper layers 
of the disk. Such settling should have little impact on 
the longest wavelength water transitions studied here, 
or on the strongest, optically thick lines traced by the 
Spitzer IRS. By creating a larger column of gas above 
the Tdust = 1 surface, the absolute fractional abundance 
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Figure 1. The Spitzer IRS high resolution (R=600) spectrum of TW Hya from 10-35 /im; the trace at bottom presents the data after the 
subtraction of a spline contin uum fit (and scaled by a factor of 3 to show the emission lines more clearly) . No water emission is detected 
in the 10-20 p,m SH modules (Najjta et al.|[20ToT ), but the 20-35 fim LH spectrum shows clear evidence for water vapor. The short vertical 
ticks indicate water vapor rest wavelengths (HITRAN 2008 database) for lines with intensities >5xl0~ 22 cm -1 / (molecule cm -2 ). 
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Figure 2. Left panel: A 2-D model fit to the TW Hya SED using RADMC and the parameters in Table[T] The stellar spectrum ( p eakin g 
near 1 fi m) is gene r ated f rom a Kurucz model w i th T e ff = 4000 K. The photometric points arc taken from Rucinski & Krauttcr (1983), 
2MASS. [Low et al.l H2005I ). and IWcintraub et alj H1989D . while the 5-35 /im Spitzer IRS spectrum is that displayed in Figured] Middle 
panel: Fiducial gas temperature profile. Right panel: Fiducial gas density profile, with the A v =l height depicted by the solid line. 



of water derived by a well mixed LTE model likely pro- 
vides an upper bound, but the inferred radial structure 
in the water vapor column density should be fairly robust 
against changes to the dust distribution. 

Above a certain altitude in the disk, where the en- 
vironment becomes exposed to the ambient radiation 
field, the gas can also be thermally decoupled from the 
dust. The gas will adopt a temperature profile that is 
a balance among heating processes, such as mechanical 
heating via accretion or that driven by photodissocia- 
tion or the photoelectric electric effect, and cooling rates 
driven by atomic, molecular, and dust grain em ission 
(|Glassgold et al.ll200i iKamp fc Dullemondl 120041) . Gen- 
eral gas temperature profiles can be estimated using de- 
tailed thermo-che mical models (e.g. iWoitke et al.1 120091 : 
iNaiita et al.ll20lH) . However, such profiles can be highly 
dependent on input assumptions, and interdependencies 
between model parameters and observables can be ob- 
scure. Here, we simplify the process by deriving the verti- 
cal gas temperature structure in the inner disks using the 



rotational ladder from M-band CO P-branch v = 1 — 
emission lines. Due to its (photo) chemical robustness, 
the abundance of CO is predicted to be fairly constant 
at nco/riHo ~ 1.2 x 10 -4 in regions warmer than 20 K 
(|Aikawa et al.l fl 996). A high CO abundance is expected 
to persist even for a depleted inner disk, su ch as that of 
TW H ya. Indeed, in the chemical model of INaiita et al.l 
(|2011|) . the abundance of CO rises to ~ 10 -4 once the 
vertical column density of H2 reaches 10 21 cm -2 for radii 
beyond 0.25 AU (physical densities are >10 9 cm -3 ). 

Since the heating of the inner disk gas is driven by 
X-ray and FUV photons from the central star and stel- 
lar ac cretion flow (|Kamp fc Dullemond 2004; Gorti et all 
120111) , we assume that gas and dust temperature become 
decoupled in regions where the radial optical depth for 
visible photons is less than unity. That is: 
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Figure 3. Model CO rovibrational fluxes from decoupled gas/dust 
temperature models compared with the observational da ta (de- 
picted as crosses with error bars, from ISalvk et alj|2007ft . The 
models show three different values of <5T, the temperature differ- 
ence between gas and dust, with values of 120 (diamonds), 135 
(squares) and 150 K (circles). 

where is dust temperature in the disk at (r, 6) in 
spherical coordinates, A v is the visual extinction along 
the radial direction 8, and 8T is a free parameter that 
describes the gas-dust temperature difference, to be de- 
termined through fits to observational data. 

The v=l -0 P(f-f2) CO lin e fluxes at 4.7 /Ltm are 
taken from ISalvk et al.l (|2007|) . while the model spec- 
tra are generated by t he raytracing code RADLitc 
(jPontoppidan et al~ll2009f) . The line s are ~7.5km/s wide 
(FWHM) (jPontoppidan et al.ll2008h . consistent with the 
nearly face-on orientation of TW Hya. The inner edge 
of the CO emission in the inner disk is set at = 
0.11 AU, based on SA imaging of the 4.7 /im line emis- 
sion (jPontoppidan et ahl 120081 ). The outer edge of the 
CO-emitting zone is not well constrained, and is set to 
4 AU, the dust transition radius. The model is relatively 
insensitive to the size of the outer edge, since the major- 
ity of the emission is produced at radii <C 4 AU. 

The best fit has £T=135± 15K (Figure 13. F rom M- 
band rotation diagram fits, ISalvk et al.1 ()2007l ) derive 
Nco = 2.2 x fO 18 cm~ 2 , corresponding to a vertical 
total gas column density of NH-T0 22 cm -2 . The f35K 
gas-to-dust temperature difference is consis tent with the 
detail ed gas temperature balance model of iNaiita et al.l 
(|20111 ). whose treatment yields <5T~100 K at r=0.25AU 
for an integrated column density of Nh = 10 22 cm~ 2 . 
RADLite simulations predict a line width (prior to in- 
strument profile convolution) of FWHM^6.8-9.6km/s 
for Uturb ~ 0.05 UKcpicr- This is somewhat larger than 
observed 7.5km/s linewidths, but this difference may be 
accounted for by the differe nce in disk inclination rela tive 
to that (i = 4°) derived bv iPontoppidan et al.l (|2008f >. 

The best fit 5T of 135 K applies only to the inner disk 
radii probed by CO, but similar physics will decouple the 
gas and dust temperatures near the disk surface at larger 
radial distances. To model this decoupling, we adopt 
T g = Td + 5T x e _r / 50AU for the gas temperature at large 
scale heights, a parame t erizat ion that matches broadly 
the models of iThi et al.l pOlOh for TW Hya. Again, the 



gas density structure follows that of the dust. 

We stress that such gas/dust thermal decoupling 
should have only a modest impact on a principle molec- 
ular mapping result presented here, namely the signifi- 
cant drop in the water vapor column density beyond the 
snow line. For the highest excitation lines measured by 
Herschel (and Spitzer) that trace the photon-dominated, 
and thus heated, layers in the inner disk and cavity wall, 
gas densities are at least ~ 10 9 H 2 /cm 3 . Un d er suc h con- 
ditions the simulations of iMeiierink et al.l <|2009f ) show 
that the emergent fluxes of the water lines detected here 
are within factors of two-three of their LTE values. Thus, 
the water vapor column densities should be reasonably 
well determined by LTE calculations whose temperature 
distributions are constrained by the CO M-band obser- 
vations. 

The outer disk is too cold to emit in such high excita- 
tion water lines, and as discussed further in Scction[5j the 
physical density in the layer at radii beyond 5-10 

AU is only -40 7 H 2 /cm 3 , some 100-1000x less than the 
critical densities of the mid- to far-infrared water vapor 
lines examined here. Only the lowest-excitation water 
transitions measured by Herschel principally trace the 
outer, deeper disk layers. Here the M-band CO obser- 
vations are of little use, but the lowest- J pure rotational 
lines do probe the relevant disk gas. In RADLite sim- 
ulati ons of lin e s up to J = 6-5 we find, in agreement 
with iQi et al.l ()2004h . that significant CO depletion in 
the cold disk midplane is needed to match the observed 
line fluxes. Water can be expected to deplete in a simi- 
lar manner, with the bulk of the vapor emission arising 
from an intermediate layer at densities where LTE cal- 
culations offer column density estimates good t o within 
an order-of- magnitude (jHogerheiide et al.ll2011| ). 



4. RETRIEVING THE RADIAL WATER VAPOR PROFILE 

4.1. The water line model 

Given the gas temperature and density structure for 
TW Hya, we can now constrain the water vapor con- 
tent of the TW Hya disk surface as a function of radius. 
Our goal is not to provide an exacting description of the 
volatile abundances versus radius and height, but to de- 
termine what radial distribution of water vapor is most 
consistent with the available data. To this end, we con- 
struct a step function in water vapor abundance, with 
one value, A"i nner , in the inner gas-depleted disk within 
4AU, another, A r i ng , in a ring starting at 4AU and ex- 
tending to the snow line at i? SI1 ow, and an outer disk 
abundance A oute r- That is, 



Ah 2 



A, 



X 



ing) 
outer ) 



R < 4 AU 

4 AU < R < R s 



(2) 



As described below, our calculations assume LTE but do 
account for line and continuum opacity. The Ah 2 o val- 
ues derived thus reflect the water vapor column densities 
to different depths into the disk. Inside of 4 AU and 
for the longest wavelength Herschel lines sensitive to the 
outermost radii the dust opacity permits the full vertical 
extent of the disk to be sampled. For the IRS features 
and the shorter wavelength PACS lines that sample gas 
near the transition radius, significant dust opacity limits 
the fitted column densities, and hence the -Yh.,o values, 
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to the upper layers of the disk. Freeze-out and dust/line 
opacity greatly limit access to the midplane beyond 4 
AU. 

Given the wide span in upper level energies of the lines 
considered, it is generally possible to identify subset of 
lines that co-vary. Consequently, we can explore one 
model parameter at a time to derive a best-fit model, il- 
lustrated in Figure HJ As we shall see, the simplest LTE 
model, namely a constant water abundance throughout 
the disk, is inconsistent with the data in hand. In Sec- 
tion 4.2.2, we further discuss whether the data support 
a modification to this basic structure in which the drop 
in water abundance beyond the snow line instead occurs 
over some finite, measurable distance. 

To render model spectra, we set the level populations 
to LTE. While the critical densities (n crit = 10 10 - 10 12 
cnr 3 ) of the high energy mid-IR water transitions sug- 
gest that a non-LTE treatment is needed, there is cur- 
rently no robust and fully tested non-LTE framework for 
the modeling of infrared water lines. That is, a non- 
LTE calculation is also likely to be inaccurate, given 
the uncertainty of collisional rates and the complexity 
of the transition network. Further, previous work sug- 
gests that non-LTE effects may alter line fluxes by factors 
of on ly a few for the mode r ate excitation lines d etected 
here (jMeiierink et all 120091 : iBanzatti et al.ll2012t l. which 
will preserve the qualitative aspects of our treatment. In 
setting level populations to LTE, we make it easier to 
reproduce our results and to evaluate the validity of our 
retrievals using more detailed models in the future. 

We assume a subsonic turbulent velocity broaden- 
ing of Uturb ~ 0-05 fKepierj and a constant ortho-to- 
para (O/P) water ratio of 3 :1, that consistent with th e 
mid- and far- infrared lines (Pontoppida n et all l2010bh . 
The cold gas seen by HIFI has a lower O /P ratio of 
0.7 (jHogerheiide et all 1201 lh , but the difference docs 
not substantially affect the derived abundance struc- 
ture. Model spectra are typically rendered with a grid 
of 0.3km/s, and then convolved with a Gaussian instru- 
ment response function to match the observed spectra. 
The Herschel HIFI spectra are rendered with a finer ve- 
locity grid of 0.05 km/s. 

4.2. Best-fitting model 

The best-fitting radial water column density profile in 
TW Hya is shown in Figure 01 whose spectral predictions 
are compared to the Spitzer LH data in Figure [5] We 
find that the inner optically thin region of TW Hya is 
dry, with upper limits on the vertically integrated frac- 
tional water abundance of Xi nner < 10~ 6 . The water de- 
tected by Spitzer originates in a thin ring near the disk 
surface, starting at 4AU and ending at a snow line just 
beyond this at R snow ~ 4.2 AU. At larger radii, the frac- 
tional water abundance drops abruptly by several orders 
of magnitude, to A ou t cr values near 3.5 xlO -11 . 

The uniqueness of this solution is illustrated in Figure 
O which shows how variations in the water abundance 
in different regions of the disk affect different lines. This 
differential sensitivity is the foundation of our method 
for retrieving the radial column density, and thus abun- 
dance, profile, and is the subject of the following discus- 
sion. 




1.0 10.0 100.0 

R (AU) 

Figure 4. The best-fit radial abundance of water vapor in TW 
Hya, illustrated by the vertically integrated column density profile 
(solid curve). Within 4 AU, the maximum water abundance is 
depicted by arrows. The dashed curve is the total gas column 
density. At the top of the figure, the role of multi-wavelength 
observations in tracing different radii of the disk are highlighted, 
along with representative facilities of each wavelength window. 




30 32 
Wavelength (/im) 

Figure 5. Detailed comparisons of the Spitzer IRS spectrum 
(black) with RADLite LTE water models (red ). The spectrum of 
RNO90 , one of the strongest water emitters in [Pontoppida n" et al.l 
Il2010bl) . is plotted at bottom, to guide the eye. The vertical dashed 
markers show the location of OH lines not discussed in this paper. 

4.2.1. The inner disk 

The inner disk of TW Hya - radii within 4 AU - is par- 
tially cleared out, leading to integrated vertical column 
densities that are more than two orders of magnitudes 
lower than those typical for cTTs. The non-detection 
of water lines between 10 and 20 /im, in the Spitzer SH 
range, may be a reflection of this, but even though the 
column density in the inner disk is dramatically reduced 
from the extrapolation of the outer disk to smaller radii, 
the scale height is sufficiently small that the physical den- 
sities at the A v ~l surface exceed 10 9 H2/cm 3 . Thus, 
non-LTE effects arc likely to be fairly unimportant, and 
we are able to produce meaningful upper limits on the 
fractional water abundance at r <4 AU. 

We investigated the effects on the Spitzer and Herschel 
lines when varying the inner disk abundance ratios be- 
tween Ai nncr = 10~ 7 and 10~ 4 . Interestingly, we find 
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Table 2 

Observed and calculated water line fluxes 



Type 


Transition 


A( /^m) 


E up (K) 


A coefficient (s ) 


Flux(10" 19 iy/m 2 ) 


Model(10- 19 W7m 2 ) 


0-H2O 


ll39-10oio 


17.23 


2438.8603 


9.619E-1 


< 163.7 


37.6 


p-H 2 


II29-IO110 


17.36 


2432.5234 


9.493E-1 


< 617.3 


110.4 


0-H2O 


836"7o7 


23.82 


1447.5970 


6.065E-1 








P-H2O 


981-872 


23.86 


2891.7024 


3.128E+1 








0-H2O 
0-H2O 


982-871 
845"7l6 


23.86 
23.90 


2891.7022 
1615.3501 


3.128E+1 
1.034 




l4669±364 a 


3051 


0-H2O 


1166-1055 


23.93 


3082.7639 


1.612E+1 








p-H 2 


1156-1047 


23.94 


2876.1493 


9.754 








p-H 2 


853-744 


30.47 


1807.0023 


8.923 


1 






P-H2O 


761-652 


30.53 


1749.8575 


1.355E+1 




>3750±317 a 


4940 


0-H2O 


762-651 


30.53 


1749.8506 


1.355E+1 








o-tl'2yJ 


054- '43 


0U.0 1 


louo.youo 


o.Oo / 


1 


>1883±287 a 


3393 


0-H20 


634-505 


30.90 


933.7488 


3.504E-1 




o-h 2 o 


8l8"7o7 


63.32 


1070.7757 


1.730 




; 123.2 


84.8 


o-h 2 o 


7o7"6l6 


71.95 


843.5459 


1.147 


< 127.2 


77.0 


p-H 2 o 


817-808 


72.03 


1270.3910 


3.099E-1 


< 35.4 


11.6 


0-H20 


423"3l2 


78.74 


432.1913 


4.852E-1 


< 39.1 


76.3 


P-H20 


615-524 


78.93 


781.1873 


4.501E-1 


< 54.3 


33.7 


P-H20 


322-2n 


89.99 


296.8471 


3.541E-1 


< 60.9 


45.8 


p-H 2 o 


4l3"322 


144.52 


396.4126 


3.345E-2 


< 5.9 


8.5 


o-h 2 o 


2l2-loi 


179.53 


114.3874 


5.617E-2 


< 8.6 


10.7 


P-H20 


I11-O00 


269.27 


53.4366 


1.852E-2 


3.07± 0.19 


3.24 


0-H2O 


lio-loi 


539.29 


60.9686 


3.477E-3 


3.46± 0.15 


1.89 



Note. — a. Water lines are often blended in the Spitzer IRS LH spectrum due to the limited resolution 
of R~600. The fluxes are measured by integrating the spectrum over each distinctive feature/region after a 
linear continuum subtraction. 



that inner disk fractional abundances higher than 10 -6 
significantly overpredict several water lines in the PACS 
range, especially the high excitation 63.3 /iin lines, while 
the Spitzer SH lines allow somewhat larger abundances. 
That is, the Spitzer SH spectrum tends to be consistent 
with a model in which the lack of lines below 20 /xm is 
explained by the overall low gas mass in the inner disk, 
and not necessarily a drier disk. It is the addition of the 
non-detection of PACS lines that require the fractional 
water abundance to be suppressed in the inner disk. 

4.2.2. The transition region and the snow line 

Next, we consider the origin of the water lines seen at 
20-35 /im. Figure |5] shows that the Spitzer LH lines de- 
tected are uniquely sensitive to variations in the water va- 
por content around the transition region. It is clear that 
a high water vapor column density is required to fit the 
Spitzer spectra, with an emitting area close to irx 1 AU 2 . 
Our fiducial dust/gas model yields A r i n g ~ 2.5 x 10~ 4 
above the Td us t=l surface (line opacity is also signifi- 
cant). As discussed in §5.2, the required column den- 
sity/abundance can be traded off against excitation tem- 
perature to some extent, with higher excitation temper- 
atures corresponding to slightly smaller radii. 

In the context of our transitional disk model, the most 
likely location for warm water vapor lies at the transi- 
tion radius, or 4 AU, a prediction that can be tested 
by measurements of the water (or perhaps OH) emis- 
sion line profiles using high resolution thermal-infrared 
spectroscopy. Interestingly, the cavity "wall" structure 
needed to reproduce the 10-30 SED (with a z/r value 
of 0.3) has a surface area close to that derived from the 
water emission for an inclination of seven degrees, further 
reinforcing the likely radial location of the high water va- 
por abundance. 

Further, the PACS and HIFI lines are sensitive to the 
location of the snow line (the outer edge of the water- 



rich ring). Consequently, we varied the location of the 
snow lines between 4.2-6 AU. In Figure [BJ it is seen that 
placing the snow line farther out in the disk over-predicts 
primarily the PACS lines. Hence, the detection of water 
lines by Spitzer in combination with the non-detection 
of water lines in PACS waveband provides strong con- 
straints on the location of the surface snow line. 

We explored one more aspect of the water abundance 
profile around the snow line. Because the water conden- 
sation temperature is a function of pressure and because 
of the strong vertical gradients in the disk, the local wa- 
ter vapor abundance need not precisely be a step func- 
tion in radius. In particular, can the water lines between 
40-150 fim be used to constrain how rapidly the water 
abundance, as inferred from the column density struc- 
ture, drops beyond the snow line with radius? 

Specifically, we modify the abundance step function 
with an exponential drop-off beyond the snow line: 

Xh,o(-R > -Rsnow) = (A r i n g— A outer )eV R eti J -|-A" ou tcr, 

(3) 

where i? cr r is the scale length of the abundance change. 
We consider i? c ff = 0.1, 0.5, and 1AU, and find that 
models with R c g > 0.5 AU produce more flux than is 
observed by the Herschel PACS instrument (Figure [7]) . 
This supports our original assumption of a step function, 
and we conclude that the surface snow line in TW Hya 
is radially narrow; that is, it occurs over a region that is 
a small fraction of its distance to the star. 

4.2.3. The outer disk 

An extremely low vertically averaged abundance be- 
yond the snow line, X outGr =3.5xlO -11 , is required by 
LTE fits to the HIFI ground state water lines at 267 
and 537 fim. Indeed, the strongest constraint on the 
outer disk water vapor abundance is provided by the 
ground-state lines, although significant limits are also 
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Figure 6. The sensitivity of selected water line fluxes to variations in the water vapor column density distribution. The first column 
shows the model water vapor distributions. The best-fit model is the solid red curve, while the green and blue curves are models calculated 
to illustrate the sensitivity of the various lines to changes in the water abundance at different disk radii. The numbers in the first column 
indicate the fractional (per hydrogen molecule) water abundance in (a) the optically thin inner disk, (b) the ring at the transition radius 
and (c) the outer disk. The final (d) panel shows the sensitivity of the model to the radius of the snow line. The following six columns 
show model line spectra compared to the observed spectra for TW Hya. The ladder crossed line in the 30 fim column is OH, which is not 
discussed in this paper. Please see Fig. [5] for additional OH line identifications. 

Xf^o ~ 10~ 4 . Outside of 4.2 AU, the vertically aver- 
aged water vapor column density decreases dramatically 
over a radial distance less than 0.5 AU, marking the loca- 
tion of the disk's surface snow line. Each of the datasets 
included in this study uniquely constrain one aspect of 
the distribution profile, and so we demonstrate the im- 
portance of the use of multi-wavelength datasets. Here 
we discuss some implications of our results. 



provided by slightly higher-lying transitions in the PACS 
range, in particular the 179.5 /im transition. The low 
water vapor co ntent of the outer disk is consistent with 
the analysis of iHogerheiide et"ail (|2011[ ) and the HIFI 
non-detect ion of water lines in another transitional disk, 
DM Tau (|Bergin et all 120101 ) . Specifically, our X outer 
is comparable to the global disk-averaged value derived 
bv IHogerheiide et all l(201 11) - 7.3 xlO 21 g of water in a 



1.9xl0~ 2 Mq disk, or X H2 o=2.2xlO- n - but as these 
authors stress there is likely significant vertical structure 
in the water vapor abundance and a large ice reservoir 
beyond 4-5 AU. 

5. DISCUSSION 

From an analysis of the water emission lines from the 
TW Hya disk over the 10 — 567 /im interval, both warm 
(^220 K) and cold water vapor are found to be present. 
The warm water emission most likely originates in a nar- 
row ring region between 4-4.2 AU where abundant water 
vapor carries much of the cosmically available oxygen, 



5.1. The origin of the surface water vapor 

How does this study inform the origin of the observed 
water vapor, and of water in protoplanetary disks in gen- 
eral? There are two potential sources of water: One is in 
situ gas-phase formation, while the other is grain surface 
formation in a cold reservoir - which could be the disk 
itself or the primordial protostellar cloud. In the second 
case, a mechanism is needed to transport the ice inwards 
and upwards to a location where it can either thermally 
evaporate or be photodeso rbed and observed in t he form 
of vapor (see, for example. iSupulver fc Linll2000l ). 
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Figure 7. Constraints on the sharpness of the snow-line in TW Hya. The left panel shows model water vapor abundance decreases from 
x if 2 ~ 2 - 5 x 10 ~ 4 to 3 - 5 x 10 -11 with three different radial scale lengths: 0.1 AU (dot line), 0.5 AU (solid line) and 1 AU (dash line). The 
right panel shows the modeled LTE line fluxes over-plotted on Herschel PACS spectra. The line styles are the same as in the left panel. 



In the warm inner disk (r< lAU, z/r<0.2), H2O 
can be formed in the g as phase (jGlassg old et all [20091 : 
iBethell fe Berghl [2009I) . even in the presence of sub- 
stantial photolyzing radiation. The key is a sufficient 
supply of H2, which drives the formation of water via 
successive hydrogenation of atomic oxygen in reactions 
with significant activation barriers. More precisely, we 
find that using the photodissociation rate calculated only 
at Lya (which dominates the 3xlO~ 3 L^ FUV contin- 
uum flux), and water reaction rates from lBaulch. D. HI 
(fl97l - k O H=3.0xl0- 14 Te- 4480 / T cm 3 molecule"^" 1 
and k H2 o=3.6xI0- 11 Te- 2590 / T cm 3 molecule^s" 1 
this simplified chemical model predicts water abundances 
near 10" 4 when T g >200K, n 9 >10 9 cm" 3 . These condi- 
tions are satisfied in the inner disk and the cavity wall, 
but not necessarily in the photon-dominated upper layers 
of the outer disk. 

Thus, production of water via gas-phase reactions can 
potentially explain the high abundance of ^220 K water 
vapor observed inside the TW Hya snow line, even if the 
initial conditions are atomic, and only until the elemental 
abund ance of oxygen is dep leted (~ 5 x 10 -4 relative 
to H, iFrisch fe Slavinl 120031). Indeed, modern thermo- 
chemical models dWoitke et all 120091: iNaiita et~aLll20lH 
iWalsh et al.ll2012l ) predict both that much of the oxygen 
is locked up in water vapor in the warm inner disk and 
at large scale heights where the gas is heated by UV and 
X-ray photons. 

The best fit water content just inside the TW Hya 
snowline is not sufficiently high to distinguish between in 
situ gas-phase formation and evaporating icy bodies, but 
is consistent with thermo-chemical models. Well inside 
the transition radius, the water abundance drops by two 
orders of magnitude, in conflict with this simple model. 
However, the models by Woitke et al. and Walsh et al. do 
not include the depleted inner region of TW Hya, making 
direct comparisons difficult for the innermost disk. 

The warm gas-phase chemistry will not operate in 



the outer disk. Here , the primordial surface ice chem- 
istry will dominate (iTielcns fc Hagenl 119821 ) , resulting 
in a huge reservoir of icy dust and larger bodies from 
which water molecules must be thermally evaporated 
or desorbed by high-ener g y par ticles and photons. As 
discus sed in iBergin et al.l (|2010f ) and iHogerheiide et all 
(| 2 1 lh . models of photodesorption from icy grains find 
that a depletion of icy grains from the outer disk sur- 
face, presumably by setting to the midplane, is needed 
to reproduce the low observed outer disk water vapor 
abundances. A zeroth-order expectation is that the wa- 
ter vapor abundanc e across the snow-line can be treated 
as a step function (jCicsla fc Cuz zi 2006), with high in- 
ner disk water vapor abundances created by a mixture 
of gas-phase reactions and evaporating icy bodies and a 
low outer disk abundance maintained by desorption. 

We find a very low (vertically averaged) outer disk wa- 
ter abundance of a few times 10 -11 per hydrogen, con- 
sistent with previous Herschcl-HIFI observations, but in 
conflict with the static thermo-chemical models, which 
indicate surface abundances of 10~ 7 — 10~ 9 per H. We 
interpret this as further evidence for settling of icy grains 
in the outer disk of TW Hya. 

Another potential test for icy-grain contributions to 
the observed disk water vapor content would be a mea- 
surement of the O/P ratio. For gas phase formation of 
water at >200K, an O/P ratio of 3 is expected, while the 
O /P ratio for evaporation or sputtering can be consider- 
ably less than 3 for grain surface temperatures less than 
~ 40 K. However, for higher temperature ices formed in 
regions close to the snow line, the O/P ratio may still 
be near 3. Thus, detailed measurements of the O/P 
ratio will be needed to trace the origin of the snow- 
line water vapor, such as could be obtained, for exam- 
ple, with next generation ground-based (VLT-VISIR) or 
space-based (JWST-MIRI) thermal-IR spectrographs. 
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Figure 8. The x 2 surface for LTE slab model water fits for TW 
Hya as a function of N and T, with an effective radius R=1.0 AU. 
The contour plot is based on Spitzer IRS LH spectra, with best-fit 
parameters of T = 220 K and A r =10 18 cm -2 . The shadowed area 
depicts the parameter space that is excluded by the water upper 
flux limits from the Spitzer IRS SH spectrum. Empty squares 
depict the best LTE slab model fit results for cTTs with wate r 
detections in the Spitzer-IKS wavelength range JSalvk et al 1 120111) . 

5.2. The water vapor abundance distribution in 
transitional disks 

The Spitzer IRS spectrum reveals that the water vapor 
emission from TW Hya is different from that of typical, 
but less evolved, protoplanetary disks around solar-type 
stars. A comparison of the basic properties of the TW 
Hya water emission is shown in Figure |8l where an LTE 
slab model (disk-averaged) temperature and column den- 
sity a re compared to those of cTTs disks (jSalvk et al.l 
1 2 1 lh . This is consistent with the detailed radial abun- 
dance structure derived for TW Hya. 

Although the relative contributions to the observed 
water vapor in TW Hya from gas-phase reactions and 
evaporating icy planetesimals is unknown, it may be the 
ideal target to observe the evaporation process as its rel- 
atively small 'hole' size allows the snow line to reside in 
the optically thick outer disk. In contrast, many of the 
known t ransitional disks have transition radii of 10 AU or 
greater ([Andrews et al J 120111 ). and the snow- line would 
be located well inside the optically thin region of the disk. 
They also typically lie two or three times more distant 
than TW Hya. Therefore, comparisons with additional 
transitional disks may need to wait for the availability of 
more sensitive instruments. 

In contrast t o the ~1 AU snow line found by 
IMeiierink et~aTl ((2009T) for cTTs disks, the TW Hya 
snow-line is located at a significantly larger radius, in 
spite of its lower luminosity as compared to DR Tau and 
AS 205 N. In the case of transitional disks, however, the 
exposure of the transition radius to nearly unobstructed 
stellar light and the 'face on' nature of the cavity wall 
at the transition radius should heat the disk more read- 
ily at large radii and push the snow-line outwards. One 
interesting consequence of this property is that if giant 
planets form and clear out a hole in the disk, creating a 
transitional disk, the subsequent migration of the snow 
line to larger radii could slow the growth of new plane- 
tary cores over a larger range of radii, out to the location 



of the new snow line. 

At the same time, new planctcsimal formation may 
be catalyzed near the new snow line due to enrichment 
of water by a combination t he radial cold finger effect 
of lStevenson fc Luninel ( 1988) and inwar ds radial migra- 
tion of icy bodies ( Ciesla fc Cuzzil [20061 ). The observed 
high abundance of water around 4 AU in TW Hya is 
supportive of a scenario in which new icy planetesimal 
formation, and perhaps even giant planet formation, is 
ongoing at this location. 

Another interesting feature of the distribution of water 
vapor abundance in TW Hya is the relatively dry inner 
disk. The water vapor abundance in the innermost disk 
(< 0.5 AU) is actually unconstrained by the Spitzer data 
due to beam dilution, but for r >0.5 AU the abundance 
upper limit is two orders of magnitude below that esti- 
mated in cTTs disks, and the lack of water in inner disk 
can be explained by two possible scenarios. 

First, the inner disk really is empty between 0.5-4 AU, 
a scenario co n sisten t with near-IR interferometric data. 
lAkeson et alj (|2011[ ) found an inner optically thick ring 
around 0.5 AU plus an outer optically thick disk starting 
at 4 AU can match the full suite of near-IR interferomet- 
ric data and the general SED shape of TW Hya. Recent 8 
- 18 /im speckle imaging, however, has suggested a more 
continuous dust distribution out to 4 AU and p erhaps 
the presence of a companion (| Arnold et al.|[2012T ). 

Further, analysis of the ALMA CO J=2-l and 3- 
2 Science Verification observations of TW Hya has re- 
yealed gas signatures down to radii as close as 2 AU 
(|Rosenfeld et all [2012] ). Thus, the more likely scenario 
for the structure of the TW Hya disk is that outlined 
in Figure [4] namely an inner disk in which CO remains 
abundant since the kinetic temperature is much higher 
than that needed for freeze-out (~ 20 K). Water in the 
inner disk can still be subjected to some depletion, or 
in the case that a companion gates the accretion flow 
into the inner disk any water-ice rich grains that settle 
to the (outer disk) mid-plane may be blocked from fur- 
ther inward migration. If this is the case, there should 
be a large difference in the gas phase C/O ratios be- 
tween the inner and outer parts of the disk. Thus, the 
quantity of gas along with its C/O content are poten- 
tially probes that can distinguish whether a given tran- 
sition disk has been sculpted by planet formation or other 
mechanism s (grain growth or photo-evaporation, for ex- 
ample; see lNajita et al.ll2011l for a further discussion of 
how changes in the C/O ratio can impact disk chem- 
istry). 

A final aspect of these observations is the short scale 
length (<0.5AU) over which the dramatic water vapor 
decrease occurs, a distance in conflict with scenarios that 
consider only freeze-out. Such a radial profile may be 
further evidence for a vertical cold finger effect, in which 
the surface snow-line location reflects that at the m id- 
plane, perhaps due to mixing (|Meiierink et al.|[2009l) . 

5.3. Molecular mapping with multi-wavelength spectra 

We have demonstrated here how multi-wavelength 
spectroscopic data can be used to reconstruct molecu- 
lar abundances, even when the observations are spatially 
(and often spectrally) unresolved. The snow line differ- 
ences inferred for the tiny sample observed to date are in 
accord with the changing physical conditions that pas- 
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sive cTTs disks experience as they evolve. Much more 
detailed constraints and comparisons to models will be- 
come available over the coming years by exploiting the 
significant infrared database that exis ts for stars+disks of 
varying mass and evolut ionary state (jPont oppida n et al.1 
l2010bHSalvk et al.ll201l| ) along with new observations at 
longer wavelengths. 

Because protoplanetary disks are complicated struc- 
tures, the multi-wavelength molecular mapping method 
should be used with some caution, however. The flux 
of each line will arise from a range of disk regions, with 
lines at various wavelengths characterizing different radii 
and/or vertical depths, and an overview of several of the 
key water vapor tracers for the best fitting TW Hya 
model is presented in Fig. [9l As demonstrated in this 
work, it can take considerable effort to build realistic 
models for an individual source, and efforts herein were 
aided by a large set of ancillary observations of this well- 
studied disk. Nevertheless, it is important to stress that 
relatively fewer uncertainties are introduced in this type 
of modeling than with the use of full thermal-chemical 
disk models, in which chemical abundances are sensitive 
to a wide array of model parameters, including FUV and 
X-ray fluxes, accretion rate, disk viscosities and trans- 
port rates, dust-versus-gas settling ge ometries, grain- 
surfa ce reaction rates, and so on (e.g. iHcinzellc r et al.l 
1201 If) . Further, studies of large disk samples will require 
fast, robust models. We therefore advocate for the direct 
measurement of chemical abundances with relatively few 
free parameters, such as is described in this work, to be 
used in conjunction with physical intuition derived from 
the more complex thermal-chemical models. 

It is interesting to note that the precision of derived 
abundance profiles is primarily determined by the avail- 
ability of high-quality spectra across a wide range of ex- 
citation energies, as well as by the knowledge of the over- 
all disk structure, especially the gas temperature. Thus, 
more precise measurements can be obtained as spectro- 
graphs with higher spectral resolution and sensitivity be- 
come available, irrespective of the instrumental spatial 
resolution (provided the disk structure is reasonably well 
characterized) - though of course instruments that pro- 
vide detailed kinematic profiles of lines with a significant 
range of excitation energies will provide the most strin- 
gent probes. This technique thus provides an ideal means 
to study chemical structure at the size scales relevant to 
planet formation. 

With a properly chosen spectral suite, we have demon- 
strated that it is possible to probe molecular abundance 
distributions on AU scales, and it is worth emphasiz- 
ing that this method is highly complementary to the 
capabilities of ALMA. At its longest baselines (16km), 
the full ALMA will be able to resolve a nearby disk 
(140 pc) on AU scales at its shortest operational wave- 
length, ^400 /Um, in dust emission. However, due to the 
quantum-limited nature of heterodyne receivers and the 
available collecting area, it will be a severe challenge even 
for ALMA to robustly image the warm molecular gas in- 
side of 10 AU - the principle formation region of terres- 
trial and giant planets according to core accretion theory. 
Thus, for the foreseeable future, multi- wavelength meth- 
ods (which can incorporate ALMA spectral data cubes) 
offer perhaps the best means of deriving the molecular 
abundance patterns in planet-forming environments. Ta- 
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Volatile line tracers at different disk radii 
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NIRSPEC, CRIRES 
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1-100 


Herschel PACS 


200- 30000 
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HIFI, ALMA 




1 10 1 00 

R (AU) 



Figure 9. The radial and vertical locations bounding water vapor 
emission from the TW Hya disk. The boxes mark the 15% and 85% 
cumulative radial line flux limits (vertical lines) and the heights 
where 15% and 85% of the line flux arises from each vertical column 
(horizontal lines). The gas density of the disk is shown in greyscale, 
and the location of the gas/dust decoupling transition (T g = 

T,j) is shown by the dotted line. As a guide to the eye, the uh 2 = 
10 9 cm -3 density contours are labeled by the solid line. Excitation 
and dust/line opacity are principally responsible for the lack of 
sensitivity to the cold, dense disk midplanc (c.f. Figure[2]l. 

ble [3] presents the general selectivity of each wavelength 
window for disks around Sun-like stars. 

Although current observational S /N levels limit chemi- 
cal abundance studies to a few species, this technique can 
in principle be applied to the study of any molecule whose 
transitions cover sufficient excitation space. The gas 
phase distribution of condensible species provides partic- 
ularly valuable information on the temperature, density, 
chemical and transport patterns in disks; and used to- 
gether, observations of a wide variety of molecules are 
best suited to the isolation and thus understanding of 
particular disk processes. To combine both infrared and 
(sub)mm data on the inner/outer disk, polar species with 
greatly differing sublimation temperatures are preferred. 
The most obvious candidates, beyond water vapor and 
OH, are CO and HCN, which are widely detected in sur- 
veys of protoplanetary disks. The former provides a fidu- 
cial chemical marker that only becomes significantly af- 
fected by freeze-out at very low temperatures, while the 
abundance drop for HCN across the snow line may well 
be sensitive to the C/O ratio in the gas. 

6. CONCLUSIONS 

This paper uses multi-wavelength spectra to probe the 
water vapor distribution in protoplanetary disks, and 
demonstrates the first application of this method to in- 
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vestigate the water vapor abundance from 0.5-200 AU in 
the transitional disk TW Hya. Our modeling shows that 
there is a narrow region between 4-4.2 AU where water 
vapor is warm (^220 K) and optically thick, resulting in 
a high abundance (X# 2 o^l0~ 4 ). Outside the snow-line, 
the water vapor column density in the disk atmosphere 
decreases dramatically over a scale length of less than 
0.5 AU due to freeze out, resulting in a vertically inte- 
grated vapor abundance of X^o^lO -11 at large radii. 

Both current and near-term astronomical facilities 
from the near-IR to radio frequencies can generate a rich 
molecular data set of the emission lines from protoplan- 
etary disks. We expect the multi-wavelength molecular 
mapping method will soon be applied to a substantial 
ensemble of disks and molecular species. 
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